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Abstract 

One of the main channels which allows for a large rate of neutralino dark matter 
annihilation in the early Universe is via the pseudoscalar Higgs A-resonance. In this 
case, the measured dark matter abundance can be obtained in the minimal supergravity 
(mSUGRA) model when tan/3 ~ 50 and 2m~^ ~ m^. We investigate the reaction pp — > 
b<fi — > bfj, + pT + X (where (j> = A or H) at the CERN LHC where requiring the tag 
of a single 6-jet allows for amplification of the signal-to-background ratio. The rare but 
observable Higgs decay to muon pairs allows for a precise measurement of the Higgs boson 
mass and decay width. We evaluate signal and background using CalcHEP, with muon 
energy smearing according to the CMS detector. We find that the Higgs width (r^) can 
typically be determined with the accuracy up to ~ 8% (~ 17%) for ~ 400 (600) 
GeV assuming 10 3 fb _1 of integrated luminosity. Therefore, the pp — > bcfi — > b^ pT + X 
process provides a unique possibility for direct Ta measurement at the LHC. While the 
Higgs width is correlated with the parameter tan/3 for a given value of m^, extracting 
tan/3 is complicated by an overlap of the A and H peaks, radiative corrections to the b 
and r Yukawa couplings, and the possibility that SUSY decay modes of the Higgs may 
be open. In the case where a dilepton mass edge from Z2 — > Z\ is visible, it should 
be possible to test the relation that 2m -x ~ m^. 
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1 Introduction 



The lightest neutralino Z\ of i?-parity conserving supersymmetric (SUSY) models is often 
touted as an excellent WIMP candidate for cold dark matter (CDM) in the universe pQ. How- 
ever, in SUSY models where the Z\ is mainly bino-like, the natural value of the relic density 
h 2 is in the 1-100 range [2], which is far beyond the WMAP observation [3], 

n CDM h 2 = Pcdm/Pc = tt CDM h 2 = 0.1123 ± 0.0035 68% CL, (1) 

where h = 0.74 ± 0.03 is the scaled Hubble constant jl]. To gain accord between theory and 
observation, special neutralino annihilation mechanisms must be invoked. These include: (i). 
co- annihilation (usually involving Z\ with a stau [5], stop j6] or chargino [7]), (ii). temper- 
ing the neutralino composition [8] so it is a mixed bino-higgsino (as occurs in the hyperbolic 
branch/focus point (HB/FP) region of mSUGRA [9]) or mixed bino-wino state or (Hi), anni- 
hilation through the light (h) or heavy Higgs boson resonance (A and/ or H) [TU] . 

In this paper, we are concerned with testing the latter annihilation mechanism, which occurs 
if 2m^ ~ rriA- The A-resonance annihilation mechanism already occurs in the paradigm 
minimal supergravity (mSUGRA or CMSSM) model, which serves as a template for many 
investigations into SUSY phenomenology. The mSUGRA parameters at the grand unification 
(GUT) scale include 

mo, m±/2, A , tan/3 and sign(fi), (2) 

where ttlq is a common scalar mass, mx/2 is a common gaugino mass, Aq is a common trilinear 
term and tan/3 is the ratio of Higgs field vacuum expectation values (VEVs). The superpoten- 
tial Higgs mass term \x has its magnitude, but not sign, determined by radiative breaking of 
electroweak symmetry (REWSB), which is seeded by the large top quark Yukawa coupling. 

In mSUGRA, as tan /3 increases, the b- and r- Yukawa couplings - /& and f T - also increase, 
and in fact their GUT scale values may become comparable to ft for tan (3 ~ 50. In this case, 
the up and down Higgs soft masses m 2 Hu and m 2 H run under renormalization group evolution to 
nearly similar values at the weak scale. Since at the weak scale m\ ~ m 2 H — m 2 Hu [11], we find 
that as tan /3 increases, the value of tua decreases [12], until finally the condition 2m ^ ~ is 
reached, whereupon neutralino annihilation through the A-resonance may take place. Another 
condition that occurs at large tan /3 is that since the b- and r Yukawa couplings are growing 
large, the partial widths T(A — > bb) and T(A — > rf) also grow, and the A width becomes 
very large (typically into the tens of GeV range). In this case, a wide range of parameter 
space actually accommodates Z\Z\ annihilation through A, H, and the value of 2m^ may be 

a few partial widths off resonance since in the relic density calculation the Z\Z\ annihilation 
rate times relative velocity must be thermally averaged. The question we wish to address here 
is: how well may one identify the cosmological scenario of neutralino annihilation through the 
heavy Higgs resonance via measurements at the CERN LHC? 

Since the 6-quark Yukawa coupling increases with tan/3, so do the Yukawa-induced Higgs 
production cross sections such as bb — >■ A, bg — > bA and gg, qq — > bbA. The presence of 
additional high pt 6-jets in the final state for the second and third of these reactions allows 
one to tag the fe-quark related production mechanisms, and also allows for a cut which rejects 
SM backgrounds - which don't involve the enhanced b Yukawa coupling - at low cost to signal. 
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The second of these reactions, which is tagged by a single fe-jet in the final state, occurs at an 
order of magnitude greater cross section than bbA production at the LHC [T5] . 

The A and H Higgs bosons are expected to dominantly decay to bb and rf final states. 
Then, the bbb or brf modes offer a substantial LHC reach for A and H, especially at large 
tan/3 [TJ]. Along with these decay modes, the decay A, H — >• has been found to be 

very useful [15] . Since the Yukawa coupling constant also increases with tan/3, this mode 
maintains its branching fraction - typically at the 10 -4 level - even in the face of increasing 
A — > bb partial width. It also offers the advantages in that the two high p? isolated muons are 
easy to tag, and the reconstruction of the invariant mass of muon pair, m^^-, allows a high 
precision measurement of the A mass and width, Fa- In fact, is was shown in Ref. [16] that the 
LHC discovery potential for A — > /i + /i~ is greatest in the b[i + fi~ mode at large tan /3, compared 
to fj, + fi~, or bbfi + fi~ . We will adopt the pp — > bA, bH + X production mode along with decay 
to muon pairs as a key to explore neutralino annihilation via the A-resonance in this paper. 

In Fig. [1], we show the leading order cross section for pp — > b<p — > bfi + fi~ + X production 
versus tua at LHC with -y/i = 14 TeV. We show curves for <fi = A and H, and for tan — 10 
and 55. Several features are worth noting. 

• The cross sections for A and H production are nearly identical, except at very low tua 
values, where substantial mixing between h and H occurs. 

• The total production cross section increases by a factor of ~ 40 in moving from tan — 10 
to tan/3 = 55. This reflects the corresponding increase in 6-quark Yukawa coupling /&, 
and goes as f£ in the total production cross section. 

• In spite of the small A, H — >• fi + fi~ branching fraction of ~ 10 -4 , the cross section for 
b/j, + n~ production via the Higgs remains large, varying between over 10 2 fb for low 

to ~ 10 _1 fb for rriA ~ 1 TeV when tan/3 is large. For LHC integrated luminosities (L) 
of order 10 2 — 10 3 fb -1 , these rates should be sufficient at least to extract the A and/or 
H mass bump. 

• In addition, the factorization scale and the renormalization scale are chosen to be [if = 
fin = m^/A with = A,H. This choice of scale effectively reproduces the effects of 
next-to-leading order (NLO) corrections [13] . 

The remainder of this paper is organized as follows. In Sec. EJ we review neutralino annihi- 
lation via the A-resonance in the mSUGRA model, and the reach of LHC for A, H fi + fi~ in 
mSUGRA parameter space for various values of integrated luminosities. In Sec. [3l we present 
our methods and results from Monte Carlo simulations for A, H production and decay to 
muons. In Section 4, we present our strategy to extract Higgs masses (m^ H ) and Higgs widths 
(Ta,h), and show the expected precision that LHC might be expected to attain in measuring 
m>A and Fa- In Sec. [5] (Conclusions), we comment on how these measurements will help as- 
certain when v4-resonance annihilation might be the major annihilation reaction for neutralino 
dark matter in the early universe. 
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Figure 1: The total production cross section for pp — > b(p — > bjj, + n + X versus in fb at 
LHC with s/s = 14 Te V. We show results for = A and H , and for tan/3 = 10 and 55. 



2 The A-resonance annihilation region in mSUGRA 

In this section, we would like to map out the portions of the 74-resonance annihilation parameter 
space which are potentially accessible to LHC searches. Figure [2] shows our results in the 
(mo, mi/2) plane of the mSUGRA model for A = 0, tan/3 = 55 and ji > 0. The green-shaded 
region has a relic densit}0 of 0.1 < Q^Ji 2 < 0.12, while the yellow-shaded region has ^^_h 2 < 
0.1. The red-shaded region has too large a thermal neutralino abundance ^^^ 2 > 0.12, and so 
is excluded under the assumption of a standard cosmology with neutralino dark matter. The 
gray region is excluded because either REWSB breaks down (right side), or we find a stau 
as the lightest SUSY particle (LSP) (left side). The blue shaded region is excluded by LEP2 
searches for chargino pair production, i.e. < 103.5 GeV. 

The A-resonance annihilation region is plainly visible on the plot. We also show the SUSY 
reach of the CERN LHC assuming yfs = 14 TeV and 100 fb -1 of integrated luminosity, taken 
from Fig. 5 of Ref. [T9] . The LHC reach is mainly determined by the total cross section for gg, 
gq and qq production, followed by their subsequent cascade decays [20] into final states with 
multi-jets plus multi-isolated leptons plus missing transverse energy (MET). A hypercube of 
cuts is examined to extract signal and background rates over a variety of cascade decay signal 
channels. We see that with L = 100 fb -1 , LHC can nearly cover the entire A-funnel. Doubling 
the integrated luminosity would allow for complete exploration of this DM-allowed region. 

1 The neutralino relic density is computed with the IsaReD [T7] subroutine of Isajet 7.80 [T5] , 
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Figure 2: The vs. m\/2 frame of the mSUGRA model for tan/3 = 55, Aq = 0, fi > 
and m t = 173.3 GeV. The green and yellow regions provide a thermal neutralino abundance in 
accord with WMAP measurements of the dark matter density. We also show contours of tua 
and r 'a/tha, arid also show the LHC reach for SUSY with 100 fb^ 1 and \/s = 14 TeV. 



Meanwhile, much of the HB/FP region is inaccessible to LHC searches, although it should be 
completely covered by future WIMP searches by Xe-100 and Xe-l-ton experiments [2"T||2"2"]. 
We also show the contour where > m e - s and where < m-^ + Mz- In the former 

region, the decay Z 2 — > ee^ — > e + e~Zi will be kinematically open while in the latter region 
the 3-body decay Z 2 — > Zie + e~ should be visible. In either case, the dilepton mass edge m e + e - 
should provide information on m^r and [23] . 

Next, we would like to know how much of the A- funnel region is open to heavy Higgs 
detection in the A — > mode. A parton level study has been performed in Ref. [16] for 

rrtA values up to 600 GeV. Here, we wish to extend these results to much higher rriA values. In 
Ref. [IE], the maximal reach for the A, H — >■ mode was found to be in the pp — > b<p + X 

channel, where <fi = A or H. 

The study in Ref. [16] evaluated pp — > b(f> — >■ + X production against SM backgrounds 

coming from bg — > bfifi, from gg, qq — >■ bbW + W~ and from gb — > bW + W~ (followed by 
W — > fiu^ decay). They required the presence of 

• two isolated opposite-sign muons with p T > 20 GeV and pseudorapidity < 2.5 , 
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• one tagged 6-jet, with p T (b — jet) > 15 (30) GeV and \r] b \ < 2.5 and 6-jet detection 
efficiency of e b = 60 (50)% for low (high) integrated luminosity regimes, 

• MET < 20 GeV (40 GeV), to reduce backgrounds from ti production for low (high) 
integrated luminosity regimes. 

The number of signal and background events within the mass range ± AM^+^- was exam- 
ined. Here, AM^+^- = 1.64 [(r^/2.36) 2 + cr^] 1 ^ 2 , with equal to the total width of the Higgs 
boson, and a m was the muon mass resolution, taken to be 2% of m^. The signal is considered to 
be observable if the lower limit on the signal plus background is larger than the corresponding 
upper limit on the background with statistical fluctuations 



L(a s + a b ) - N^L{a s + a b ) > La b + Ny/La b (3) 

or equivalently, 



N 



N + 2JLa b 



(4) 



Here L is the integrated luminosity, a s is the cross section of the signal, and a b is the background 
cross section. The parameter N specifies the level or probability of discovery, which is taken to 
be N = 2.5 for a 5<r signal. For a b ^> a s , this requirement becomes similar to 



N s La s 



where N s is the signal number of events, N b is the background number of events, and N$s = 
the statistical significance, which is commonly used in the literature. 

Here, at the first stage of our analysis, we repeat this calculation, although we extend the 
results to much higher values of and higher integrated luminosities. We also have evaluated 
additional possible backgrounds to make sure that their contributions are either important or 
negligible. In our analysis, we use the CTEQ6L set for PDFs [21] and the QCD scale is set equal 
to myi/4 for signal and s for backgrounds. The following backgrounds have been evaluated with 
the respective i^-factors applied to take into account Higher order corrections: 

• 99 + QQ ~^ W + W~bb — > fi + yTvv [K = 2): this is the dominant background coming 
mainly from ti production and decay, 

• bg — > W + W~b — > fi + fi~uu (K = 1.3): this background is typically at least one order of 
magnitude below the first one, 

• bg — > bfi + n~vu (K = 1.3): this background is of the same order of magnitude as previous 
one, 

• 67 — > bn + fi~uu (K = 1.3): this background is several times lower than the previous one 
and can be considered as a subdominant one, contributing to the total background at 
the percent level. It was evaluated using the photon distribution function of the proton 
available in CalcHEP, 
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• eg — >■ cfi + fi~uu (K = 1.3): this background is of the same order as the previous one 
and therefore, again contributing to the total background at the percent level. It was 
evaluated using a mis-tagging probability for c-jet equal to 10%. 

• qg — > qfi + fi~uu (K = 1.3): also, this background is of the same order as the previous 
one, contributing to the total background at the percent level. It was evaluated using 
mis-tagging probability for light g-jet equal to 1%. 

Fig. |3] presents pp — > <p°b — > p + p~b + X signal rates versus m^, where 0° = A,H,h after 
application of kinematical cuts and efficiency of 6-tagging for tan (3 = 5 (black lines) and 
tan/3 = 55 (red lines). Results for low and high luminosity regimes are denoted by solid and 
dashed lines respectively. 




200 400 600 800 1000 

m A (GeV) 

Figure 3: pp — > cj)b — > fi + fi~b + X signal rates versus m^, where <p = A, H, h after application 
of kinematical cuts and efficiency of b-tagging for tan/3 = 5 (black lines) and tan/3 = 55 
(red lines). Results for low and high luminosity regimes are denoted by solid and dashed lines 
respectively. 

In Fig. HI we present rates for various backgrounds described above for p + fi~b signature 
versus after application of kinematical cuts and efficiency of 6-tagging for an intermediate 
value of tan /3 = 30. Results for low and high luminosity regimes are presented in left and right 



6 



frames respectively. One can see that indeed the contributions from the last three subdominant 
backgrounds discussed above are at the percent level. 




Figure 4: Rates for various backgrounds for n + n~b signature versus rriA after application of 
kinematical cuts and efficiency of b-tagging for tan (3 = 30. Results for low and high luminosity 
regimes are presented in left and right frames respectively. 

Using signal and background rates from these calculations, we derive the LHC discovery 
reach. The results are shown in Fig. |5j One should notice an important effect of the cuts 
for low and high luminosity regimes. The main effect for LHC reach comes from the MET 
cut. We require MET < 20 for low luminosity which should leave signal intact assuming that 
instrumental missing transverse momentum is under control above 20 GeV in the low luminosity 
regime. This cut significantly suppresses the leading ti background. For high luminosity regime 
we apply MET < 40 which does not affect signal but significantly increase background. The 
overall effect of high luminosity cuts is an increase of the background and decrease of the 
signal. Therefore, the discovery potential of the LHC at L = 100 fb~ is slightly lower than at 
L = 30 fb _1 . But in the region sufficiently above the border-line for LHC discovery potential 
shown in Fig. [51 say for tua = 400 GeV and tan /3 = 55, LHC at L = 100 fb _1 provides better 
statistics and significance as compared to the L = 30 fb _1 case as we show below. 

Here, we see that for L = 100 fb _1 , the reach for 60 — > bfi + fi~ at tan/3 ~ 55 extends to 
rriA — 550 GeV. For L = 300 fb _1 , the reach extends to vtla — 730 GeV, and for L = 1000 fb^ 1 , 
the reach extends to ~ 925 GeV. 

From the results of Fig. [51 we can now compare against Fig. |2] to see how much of the 
A-funnel can be explored via the A, H — >■ decay mode. To illustrate, we show contours 

of tua = 500, 750, 1000 and 1250 in Fig. |2j Thus, for 100 fb" 1 of integrated luminosity, we 
expect LHC to be sensitive to a A, H — > n + fi~ bump for about half of the A-funnel. An 
integrated luminosity of 300 fb _1 covers about three-quarters of the A-funnel, while well over 
1000 fb" 1 will be needed to cover the entire funnel region. 

In addition to measuring the value of itia,h y i a a dimuon mass bump, one may be able to 
extract information on the A, H widths from the dimuon channel, if LHC experiments have 
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Figure 5: LHC reach for the pseudoscalar Higgs A via pp — > bA — > bfi + fi~ + X in the 
rriA vs. tan/3 plane for various possible values of integrated luminosity. The 30 fb" 1 reach 
exceeds the 100 fb~ l reach because we use harder cuts in the high luminosity case. 

sufficiently good muon energy reconstruction. To illustrate the values of Fa that are expected, 
we plot contours of T^/m^ in Fig. |2j These range from about 5% for low myi ~ 250 GeV, 
corresponding to Y A ~ 15 - 20 GeV, to about 4.4% for m 1/2 ~ 1200 GeV, where Y A ~ 50 GeV. 
To better illustrate the range of Higgs widths expected in mSUGRA, we show in Fig. [H] the 
value of Ta versus m A after a scan over mSUGRA parameter space for various fixed values of 
tan/3. Here, we see that indeed as tan/3 grows, so too does V a- In fact, for a measured value of 
mA, a measurement of Ta will indicate a rather small window of allowed tan (3 values. Naively, 
one might expect a one-to-one correspondence between tan/3 and Ta for fixed m^. However, 
two effects that spread out the correlation include: (i). weak scale threshold corrections to fb 
that are large at large tan /3, and depend on the entire SUSY spectrum via loop effects [25], and 
(ii). various additional SUSY decay modes of the A and H may open up (26], depending on 
sparticle masses and mixings. For instance, if rriA > 2m^, then the decay mode A — > W^~W{~ 
opens up and contributes to the A width. Thus, models with lighter SUSY particles should 
correspond to larger Ta values for a given and tan /3 value, whereas if all non-standard 
decay modes are closed, then the lower range of Ta that is shown may be expected to occur. 
The loop corrections to fb tend to enhance fb for /i > and diminish /& for /x < 0, leading to 
somewhat separated bands for each tan (3 value. 
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Figure 6: Plot of Ta vs. vtia from a scan over mSUGRA model parameters for tan/3 = 
10, 20, 30, 40 and 50. 

3 Detailed simulations for pp — > bA, bH + X 

In this section, we present a detailed Monte Carlo study of detection of b<p — > bfi + p~ for a 
particular case study. The benchmark point we adopt is known as LCC4 in the study by 
Battaglia et ai, Ref. [27]. Some of the mSUGRA parameters and sparticle masses as generated 
by Isajet 7.80 are given in Table [TJ We use a value of m t = 175 GeV instead of 178 GeV as in 
Ref. [27J since the latest Isasugra/IsaReD code gives a relic density of ^ h 2 = 0.1 for the 175 
GeV value, and 0.16 for the 178 GeV value. We also examine later how well the value of Ta 
can be measured for benchmark point BM600 with tua = 608 GeV. 

The resolution of the dimuon invariant mass, and hence an accurate measurement of vtia 
and especially Ta, depends on the LHC detector's ability to measure the muon's momentum. 
The muon momentum is measured from its amount of bending in the magnetic field of the 
detector. Thus, for low energy muon, with a highly curved track, the muon p measurement 
should be more precise than for high energy muons, which have very little track curvature. For 
our studies, we use a CMS muon smearing subroutine, where the smearing as a function of 
\rj(p)\ is displayed in Fig. [7J for several muon p T values [2"8||2"9] . 

We begin our MC simulation by calculating bg — > bp + pT production for pp collisions at 
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325.4 


519.5 


m z, 


172.5 


274.7 


m A 


420.7 


607.9 


m H 


A OQ £ 


olz.U 


m h 


1 1 c 1 

115.1 


1171 

117.1 




or . , 1 r\ — -10 

35 x 10 


11 x 10 


BF{b 57) 


1.9 x 10- 4 


2.8 x 10^ 4 


BF(B S -> n+n~) 
Qh 2 ~ 

a(Z ± p) pb 


2.8 x 10~ 8 


1.1 x 10" 8 


0.096 


0.089 


1.1 x 10~ 8 


1.7 x 10 9 




19.1 GeV 


31.9 GeV 




19.2 GeV 


32.1 GeV 



Table 1: Masses and parameters in GeV units for Benchmark points LCC4 (with m t = 175 
GeV) and BM600 (with m t = 173.3 GeV) using Isajet 7.80. 

y/s = 14 TeV using CalcHEP [30]. The relevant Feynman diagrams are displayed in Fig. [SJ 
They include not only A and H production and decay, but also background contributions from 
7*, Z* and h production. 

In Fig. [91 we plot the invariant mass distribution of muon pairs for L = 30 fb _1 . 

For all distributions now and hereafter, we take into account detector effects of muon momenta 
resolution according to Fig. [7| using Gaussian smearing applied to the particle's momentum 
generated by CalcHEP at the parton level. What is clear from the plot is that the 7, Z — >■ 
peaks stand out; but also the A, H overlapping peak stands out well above background 

levels at m A ~ 420 GeV. 

In Fig. [101 we plot the muon distribution (solid line) and b-jet p? distribution (dashed 
line) from pp — » bA — > bfi + fi~ + X production for the LCC4 benchmark. The muon 
distribution peaks at around px ~ Tn A /2, but with substantial smearing to either side due to 
the momentum of the A. Since the fe-jets are emitted preferentially in the forward direction, the 
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Figure 7: Plot of CMS muon smearing function versus 1 77 (/x) | for various muon Pt values. 
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Figure 8: Feynman diagrams for bg — > m i/ie MSSM. 



Pr{b) distribution peaks at low values, with some smearing out to values over a hundred GeV. 
In Fig. HU we plot the muon (solid line) and 6-jet (dashed line) pseudo-rapidity distributions 
from pp — > bA — > bfi + fi~ +X for the LCC4 benchmark. The muon 77 distribution is clearly more 
central, while rj(b) is less central due to its role as an element in QCD initial state radiation. 
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Figure 9: Plot of invariant mass distribution of muon pairs from a CalcHEP MC com- 

putation using benchmark LCC4- 
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Figure 10: Plot of distribution in Pt(i^) (solid line) and pr(b) (dashed line) from a CalcHEP 
MC computation using benchmark LCC4- 
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Figure 11: Plot of distribution in t](fi) (solid line) and t](b) (dashed line) from a CalcHEP MC 
computation using benchmark LCC4- 

4 Extracting rriA,H and Va,h 

Once a dimuon mass bump has been established, the next step is to fit the invariant mass 
distribution with a curve which depends on the Higgs mass and width. A complication occurs 
because in our case the A and H masses are only separated by ~ 3 GeV, and so the two peaks 
are highly overlapping, and essentially indistinguishable. To see what this means for an ideal 
measurement, we plot in Fig. [12] for LCC4 the dimuon invariant mass from just the reaction 
pp — >■ bA — >■ bfi + fi~+X (red curve), and also the distribution from pp — > bH — > bfi + fi~+X (blue 
curve), along with the sum (black curve). A direct measurement of these idealized distributions 
of full-width-at-half-max shows indeed that T A ~ 430-410 = 20 GeV, while T H ~ 433-413 = 
20 GeV. A measure of the summed distributions provides V a,h — 433 — 410 = 23 GeV, i.e. the 
idealized width expectation expanded by the A, H mass splitting. We fit the dimuon invariant 
mass distribution from all diagrams of Fig. [8] along with muon smearing with the following 
function F of dimuon mass m and 6 fitting parameters T, M, N, a, N p i, N p2 : 



F{m; T, M, N, a, JV p i, N p2 ) = N / B(rri, T, M) x G{m! , m, a)dm! + iV p i exp(-iV p2 m), (5) 



where N is just a normalization parameter, 



B(m',T,M) 



2 r 2 M 2 



G(m' , m, a) 



it (m' 2 — M 2 ) 2 + m' 4 (r 2 /M 2 ) ' 
1 {m! — m) 2 
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Figure 12: Plot of dimuon invariant mass (M^+^-J from bA production (red), bH production 
(blue) and sum (black) for benchmark LCC4 with no smearing. 

One can see that F(m; T, M, N, a, N p i, N P 2) is a convolution of the Breit-Wigner resonance func- 
tion along with Gaussian detector smearing plus an exponentially dropping function describing 
the background shape. 

The results from the x 2 fits of signal-plus-background are presented in Fig. [TB] for different 
integrated luminosities of L =30, 100, 300 and 1000 fb _1 . The left side of the Fig. [TBI shows the 
fit to Monte Carlo data for LCC4 with pp — > bA, H — > bfi + fi~ + X production including muon 
smearing. It also shows the values of the fitted parameters (r, M, N, a, N p i, N P 2) together with 
their standard deviations according to the fit. The fit has been performed using the MINUIT 
program from CERN library which properly takes into account the correlation matrix of the fit 
parameters, which is crucial for the evaluation of the corresponding contours in the Ta vs. tua 
plane at la and 2<r confidence levels; these are shown on the right side of the Fig. [131 The black 
crosses show the width measurement assuming the Higgs masses are known to perfect accuracy. 
From Fig. [TBJ, one can see that for L =30 fb _1 of data, the statistics only provide a rough fit 
to the A, H width. On the other hand, moving to L =100 fb" 1 , our fit provides promising 
results for Ta- We see that with L =100 fb _1 , rriA can be measured to 1 GeV accuracy, or 
0.25%. Meanwhile, the A, H width is measured at T a ,h — 20 ±8 GeV, or 40% level. At higher 
integrated luminosity values of L =300 fb _1 , the accuracy on rriA,H is improved to sub-GeV 
levels and Ta,h is found to be ~ 20 ± 4 GeV, a 20% measurement. At L =1000 fb _1 , which 
might be reached in ~ 10 years of LHC running, the measurement of V a,h can be improved to 
about 20 ±1.75 GeV, or ~ 8% accuracy. The Ta,h accuracy is expected to approach ~ 7% level 
for infinite integrated luminosity, and is mainly limited by the detector muon energy resolution 
of 4%, which is actually quite close to Y AjmA- 

The results from the \ 2 fits to rriA and Ta of signal-plus-background for benchmark point 
BM600 are shown in Fig. [TT] for integrated luminosities of L =300 and 1000 fb _1 . In this 
case, the heavy Higgs masses are vtla^h — 608, 612 GeV, while the widths from CalcHEP are 
Fa,h = 31.9 (32.1) GeV. We find for 10 3 fb _1 of integrated luminosity, that Ta is extracted to 
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Figure 13: 

production 
for Monte 
smearing. 



Left : best fit of Monte Carlo data for LCC4 for pp — > bA, H — > bfi + fi~ + X 
including muon smearing. Right: corresponding contours of fit to tua and Ta values 
Carlo data for LCC4 from pp — > bA,H — > bfi + fi~ + X production including muon 
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Figure 14: Left : best fit of Monte Carlo data for BM600 for pp bA, H + X 

production including muon smearing. Right: corresponding contours of fit to tua and Ta values 
for Monte Carlo data for BM600 from pp —7 bA, H — > bfi + fi~ + X production including muon 
smearing. 



be = 28 ± 5.5 GeV, a 17% measurement. 



5 Conclusions 

In i?-parity-conserving supersymmetric models where the lightest neutralino is expected to 
be a thermal relic of the Big Bang, and to comprise the dark matter in the universe, special 
qualities are needed to enhance the dark matter annihilation rates. One compelling case is 
neutralino annihilation through the pseudoscalar Higgs resonance. Can the LHC tell us if such 
a mechanism is operative in the early universe? The crucial test here is whether the condition 
2m^ ~ rriA is fulfilled. 

A variety of techniques have been proposed for extracting the SUSY particle masses - 
including - in sparticle cascade decay events at the LHC [51] . Extracting the heavy 

Higgs masses is also possible provided that rriA is small enough and that tan (3 is large enough. 
Mass measurements of heavy Higgs decays into bb and t + t - are fraught with uncertainties 
from multi-particle production and energy loss from neutrinos. We focused instead on the 
suppressed decay A, H —7 since it allows for both highly accurate heavy Higgs mass 

and width reconstructions. Production of A and H in association with a single 6-jet offers a 
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large background rejection at small cost to signal, especially in the large tan/3 regime, where 
Higgs production in association with bs is expected to be enhanced by large Yukawa couplings. 
This is also the regime in models such as mSUGRA where neutralino annihilation through the 
heavy Higgs resonance is expected to occur. 

In this paper, we have computed regions of vtla vs. tan /3 parameter space where pp — > 
b(f) + X,(f) = H,A production followed by H, A — > should be visible for various integrated 

luminosities. We have also performed detailed Monte Carlo simulations of signal and back- 
ground for the LCC4 and BM600 benchmark points. Fits of the dimuon mass spectra allow for 
sub-percent determinations of the (nearly overlapping) H and A masses. The A, H overlapping 
widths were determined to ~ 8% (~ 17%) accuracy in the case of LCC4 (BM600) with 10 3 fb _1 
of integrated luminosity. We conclude that indeed the study of pp — > b<fi + X, <p = H, A produc- 
tion followed by if, A — > /i + [i~ offers a unique opportunity to directly measure A(H) Higgs 
width. This process also allows to measure the A (H) mass with unprecedented precision. Both 
these measurements would provide crucial information to connect the cosmological A — funnel 
scenario of dark matter annihilation with LHC data. Combining these measurements with 
SUSY particle mass measurements such as the mass edge in n%i+i- from Z 2 —> Zi£ + £~ decay 
would go a long way towards determining the parameter tan /3, and also whether or not neu- 
tralino annihilation through the A resonance (with 2m^ ~ m^) is the operative mechanism in 
the early universe to yield the measured abundance of neutralino dark matter. 
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